Aeolian sans ripples: experimental study of saturated states 
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We report an experimental investigation of aeolian sand ripples, performed both in a wind tunnel 
and on stoss slopes of dunes. Starting from a flat bed, we can identify three regimes: appearance 
of an initial wavelength, coarsening of the pattern and finally saturation of the ripples. We show 
that both initial and final wavelengths, as well as the propagative speed of the ripples, are linear 
functions of the wind velocity. Investigating the evolution of an initially corrugated bed, we exhibit 
non-linear stable solutions for a finite range of wavelengths, which demonstrates the existence of a 
saturation in amplitude. These results contradict most of the models. 

PACS numbers: 45.70.-n, 47.54.+r, 05.45.-a 



The surface of aeolian sand dunes is not smooth but is 
usually formed into regular patterns of ripples, transverse 
to the wind 1]. Their wavelength ranges from the cen- 
timeter to the meter with a constant aspect ratio (~ 4%) 
. Although many different models have been proposed 
to explain the formation and evolution of aeolian rip- 

MpUl S II II S H ES EH, few field observations 
lia | and controlled experiments 0,0,0 have been 
performed so far. By contrast with subaqueous dunes 
or ripples which result from a hydrodynamic instability 
induced by the interaction between shape and flow 
aeolian ripples are of different nature and result from a 
screening instability. When the 'saltons' - high energy 
grains - collide the bed, they eject grains of smaller en- 
ergy, 'reptons'. The windward slope of a small bump is 
submitted to more impacts than the lee slope, so that the 
flux of reptons is higher uphill than downhill and makes 
the bump amplify. 

Most of the models agree for the linear stage of the 
instability - see for a pedagogical review. They as- 
sume that the reptons remain trapped by the bed after a 
single hop of length a, distributed according to a distri- 
bution -P(a). The saltons are considered as an external 
reservoir which brings energy into the system. On this 
basis, the most unstable wavelength Ao is found to scale 
on the reptation length a — Jda a P(a). Besides, the 
most recent investigations of sand transport, both exper- 
imental [la, and theoretical [171 Il8| . indicate that 
P(a) is independent of the wind shear velocity u* , which 
implies that a scales on the grain diameter d. The ini- 
tial wavelength of the ripples Ao is thus expected to be 
independent of the wind strength. 

In the non-linear regime, aeolian ripples have been one 
of the most extensively studied systems presenting coars- 
ening. Two scenarii are evoked. The wavelength can in- 
crease by a phase instability (negative diffusion) leading 
to a continuous stretching of the pattern [J^. Alter- 
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FIG. 1: Experimental set up both in the field and in the wind 
tunnel, (a) The profile h(x) is given by the inclined laser 
sheet, (b) Autocorrelation function C(\,t) =< h(x,t)h(x + 
X,t) > — < h(x,t) > 2 of the profile. The position of the 
first peak gives the m ean wave length A. The amplitude is 
defined by A(t) = 2- v /2C(0,t). (c) Method for engraving 
periodic profiles (the lateral bars are removed once the pattern 
is engraved). 



natively, ripples may be considered as locall y in teracting 
objects that undergo successive merging (l3.l20| . In both 
cases, one predicts a perpetual increase of X(t) as the log- 
arithm of time or as a power law of small exponent. Such 
a logarithmic behaviour has been argued to be compat- 
ible with some wind tunnel data [14|. However, com- 
mon field observations clearly evidence steady patterns 
driven by a wind permanently blowing for hours. 'Fully 
developped' ripples are also reported in experiments of 
EH, but not studied in details. Some possible reasons to 
explain the discrepancies between models and field ob- 
servations are discussed in 0, and to the best of our 
knowledge, only one recent model predicts a saturation 
of ripples wavelength at long time The goal of this 
paper is to investigate experimentally the precise condi- 
tions of saturation, and to study the properties of the 
saturated states. 
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FIG. 2: Evolution of an initially flat bed. (a) Pictures at different times for field experiments u* = 1.3 Uth, (b) Spatio temporal 
diagram in the wind tunnel for u* — 1.4 u t h- Time evolution of (c) wavelength \(t) and (d) amplitude A(t) from a flat bed for 
different values of the shear velocity: u* = 1.3 u t h (A), u* = 1.4 Uth (A), u* — 1.8 Uth (O), u* = 2.3 Uth (•). Inset: time log 
scale. 



Experimental set-up - We have performed experiments 
both in the wind tunnel of the CEMAGREF in Grenoble 
(1 m wide, 0.5 m high and 4.5 m long) and in the barchan 
field extending between Tarfaya and Laayoune in Mo- 
rocco. In the wind tunnel, the air flow is well controlled 
and can reach u* ~ 0.7 m/s but the sand flux is not 
fully saturated due to the limited length of the sand bed 
(3 m). In the field, the sand supply is unlimited and 
the flux is saturated but the wind strength fluctuates. 
The grains we used in the lab - so-called Hostun sand - 
are 120 ± 40 fim in diameter and are rather angular as 
this sand comes from a quarry. The grains in the dune 
field is more round and smooth, made of quartz/lime 
mixture, and with d ~ 180 ± 35 /im. In both situa- 
tions the measured shear velocity threshold for saltation 
is Uth ~ 0.22 m/s. To measure the sand surface defor- 
mation, the same method is used both in the field and 
in the wind tunnel. A laser sheet is inclined at some low 
angle to the bed (Fig. HJi) and pictures are taken from 
above using a 2240 x 1680 digital camera. The profile 
h(x, t) is detected by a correlation method which insures 
a precision of 40 ^m. The patterns obtained are usually 
quasi-periodic and can be characterized by an average 
wavelength A and an amplitude A (Fig. \\}p). 

Instability of a flat bed - We first study the time evo- 
lution of an initially flat bed. The dynamics in the field 
is shown in Fig.|2i and the spatio-temporal diagram ob- 
tained in the wind tunnel is displayed in Fig.03. In both 
cases, the initially flat surface deforms and coarsening 
is observed. The time evolution of wavelength and am- 
plitude for experiments in the wind tunnel are plotted 
in Figs. [2t,d. Identifying a linear regime is not an easy 
task. A gentle linear instability is usually characterised 



by an exponential growth of the amplitude and a plateau 
in wavelength. Here, tiny structures with a small wave- 
length show up after a short time (~ 1 min for this run). 
Their growth in amplitude is consistent with an exponen- 
tial when A is smaller than, say, 1 mm i.e. 6 grain sizes. 
However, a peak clearly appears in the autocorrelation 
function C(A, t) only close to the end of this period. This 
tiny plateau still allows to define the initial wavelength 
Ao- As soon as they are visible, the ripples start merging 
so that the pattern exhibits coarsening. The increase of 
amplitude and wavelength progressively slows down and 
after some time (typically ~ 10 min), the pattern tends 
to 'saturate': the ripples essentially propagate without 
changing shape and amplitude anymore. We take Aoo 
as the average of A over the last two minutes. As em- 
phasized by the log scale in time, this ultimate regime 
could also be interpreted as a slow drift in A (lif. Set- 
tling this issue would require at least one more decade 
in time (few hours), which is difficult to achieve in our 
wind tunnel due to the global erosion of the bed (a sand 
layer few centimetres thick disappears after typically half 
an hour) and to the convective nature of the instability 
(the entrance conditions can influence the measurement 
zone) . 

Evolution of an initially corrugated bed - To bypass 
the problem, we have performed experiments in which 
we do not start from a flat bed but from a corrugated 
surface. A periodic ripple profile is initially engraved 
with a cardboard pattern on the bed and is then let free 
to evolve (Fig. U;). In the wind tunnel, we have been 
able to systematically study the influence of the initial 
wavelength and amplitude, the shape remaining that of 
a developed ripple. The first result concerns the case 
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FIG. 3: Evolution of an initially corrugated bed. Pictures from the field for u* — 1.3 Uth and initial wavelength (a) A = 5 cm; 
(b) A = 13.5 cm; (c) A = 30 cm. Spatio temporal diagram in the wind tunnel for it* = 1.3 Uth and (d) A = 5 cm; (e) A = 9.3 cm; 
(f) A = 19.5 cm. 



when the engraved pattern has the same wavelength Aoo 
and same amplitude as measured in 'natural condi- 
tions', i.e. starting from a flat bed. The ripples do not 
evolve at all and the pattern purely propagates (Fig.Efc), 
showing that a saturated state exists. If we start from 
the same A but with a larger or a smaller amplitude, the 
pattern converges back to the same shape (A^, X^) by 
adjusting its amplitude. If we start from the same A but 
a slightly different A, the patterns keeps this new wave- 
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FIG. 4: Stability diagram of ripples of initial wavelength A 
and amplitude A for u* — 1.3 Uth- Outside the locking 
zone (O), the pattern destabilizes by generation of harmonics 
(right) or sub-harmonics (left). Starting in the central region 
(■), A stays constant and A readapts to reach a stable steady 
solution (•). Very large amplitudes cannot be imposed as 
avalanche processes take place. Hatched region is beyond ex- 
perimental access. 



length and selects a new amplitude (Fig. 0J. We can 
thus conclude that, depending on the initial conditions, 
the pattern converges toward different stable non-linear 
solutions. The saturated state observed starting from a 
flat bed is one solution among the finite range of wave- 
lengths for which stable ripple patterns can be observed. 
If we start from a too large A, the pattern destabilizes by 
changing its wavelength, forming substructures on the 
windward slope (generation of harmonics) as shown in 
Fig. [3f- Too small A are also unstable as the pattern 
coarsens due to merging (period doubling) (Fig. |3|i) . In 
the (A, A) diagram representing the different initial con- 
ditions (Fig. , one can then determine a stable region 
where the engraved pattern does not change wavelength, 
an unstable region leading to coarsening for A < 60 mm 
and an unstable region leading to a windward slope in- 
stability. These observations in the wind tunnel hold for 
our field experiments: Figs|3^,,b,c show that the engraved 
ripples follow the same spatio temporal evolution as in 
the wind tunnel when the wavelength is smaller, close to 
or larger than the 'natural' one. Due to wind fluctua- 
tions we have not been able to systematically study the 
influence of the initial amplitude. 

Parametric study - Once we have shown that the final 
pattern observed starting from a flat bed corresponds to 
a saturated state, we can systematically study the influ- 
ence of the wind on the natural ripples characteristics. 
Fig. [SJi shows the selection of the aspect ratio for differ- 
ent wind strengths. While the aspect ratio remains con- 
stant in the field (A ~ 0.04A), the amplitude saturates 
in the wind tunnel at large u*. This discrepancy may 
be related to the unsaturated state of the sand flux. The 
initial and final wavelengths scaled by the particles diam- 
eter are plotted as a function of the wind speed in Fig.JSjD. 
By contrast with standard predictions, they both increase 
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FIG. 5: (a) A as a function of A for different wind velocities in the wind tunnel (•) and in the field (O). (b) Initial (■) 
and final (•) wavelength (wind tunnel data) and field final wavelength (O) as a function of the imposed flow velocity u». (c) 
Propagation velocity of developped ripples as a function of the wind speed for free evolving ripples in the wind tunnel (•) and 
in the field (O). Data represented by + (resp. x) come from [T^| (0), for grains of size d ~ 310 /im (150 fim). 



linearly with u* . This linear behaviour was also found by 
[l3|. Interestingly, the initial wavelength vanishes at the 
transport threshold and both Ao (it*) and Aoo(w*) are par- 
allel. The same tendency holds for field measurements, 
although more scattered. In a similar way, the rescaled 
phase velocity c 00 /\/~gd of saturated ripples increases lin- 
early with the wind (Fig.^t). Our data then suggest that 
A /d oc {u*/u t h - 1) and Coo/^/gd, cx (u*/u th - 1). In or- 
der to test the scaling with d which simply comes from 
dimensional analysis, we have also plotted in Figs. Eh, c 
data coming from and 01 The collapse is good, 
showing that d is the relevant length scale. The prefactor 
in these linear relations certainly depends on the density 
ratio Psand/Pair, which however is difficult to vary. Com- 
parison with martian ripples formed in an atmosphere 
80 times lighter [2lJ could be of interest to get deeper 
insight in the scaling. 

Conclusion - Two main conclusions can be drawn 
from this study. First, the linear dependence of the 
most unstable wavelength with the wind speed rules out 
our theoretical understanding of ripples formation and 
sand transport in general. As Ao turns out to be much 
larger than the value predicted by Anderson-like models 
(6a ~ 60d ~ 10 mm), there should exist a yet unknown 
mechanism stabilizing small wavelengths. Second, we 
have shown the existence of a whole family of stable non 
linear solutions. Future models should then be able to re- 
cover the range of wavelengths for which ripple patterns 
are stable, as well as the saturation of A when starting 
from a flat bed. 

In this work, we have studied the one dimensional dy- 
namics of the ripples along the wind direction. However, 
a complex dynamics exists in the transverse direction too 
and the formation of defects is known to play a significant 
role in coarsenin g pr ocesses and in the route towards reg- 
ular pattern [fI l2(lj . It then suggests to investigate the 
density of defects to characterize the degree of satura- 
tion of the pattern. Finally, ripples could be of practical 



interest for field experimentalist as a non intrusive mea- 
surement of the reptation sand flux. 
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